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ABSTRACT/RESUME 

NASA Goddard Space Flight Center and the European 
Space Agency (ESA) conducted a disturbance 
verification test on a flight Solar Array 3 (SA3) for the 
Hubble Space Telescope using the ESA Large Space 
Simulator (LSS) in Noordwijk, the Netherlands. The 
LSS cyclically illuminated the SA3 to simulate orbital 
temperature changes in a vacuum environment. Data 
acquisition systems measured signals from force 
transducers and accelerometers resulting from 
thermally induced vibrations of the SA3. The LSS with 
its seismic mass boundary provided an excellent 
background environment for this test. This paper 
discusses the analysis performed on the measured 
transient SA3 responses and provides a summary of the 
results. 

1 INTRODUCTION 

This year during the Hubble Space Telescope (HST) 
servicing mission 3B (SM3B), NASA astronauts will 
install Solar Array 3 (SA3) on HST. This will replace 
the nearly 7 year old Solar Array 2. The SA3 will 
provide more power with less surface area by using 
Gallium Arsenide Solar Ceils mounted on honeycomb 
panels. NASA expects the new SA3 to provide power 
to HST for the remainder of its mission. Fig. 1 
illustrates how HST will appear after SM3B. 



To qualify for flight on HST, a test was designed to 
verify that the SA3 would not disturb the HST after 
installation. NASA called this test the Disturbance 
Verification Test (DVT) of Solar Array 3. The test 
objective was to measure thermally induced vibrations 
of SA3 in a test simulated space environment. The 
events of interest were transient vibrational responses. 
It was mandatory to measure these vibrations if they 
had been large enough to cause Loss of Lock (LOL) on 
the HST guide stars. [1] 

To conduct this test, NASA Goddard Space Flight 
Center (GSFC) used the European Space Agency 
(ESA) Large Space Simulator (LSS) in Noordwijk, The 
Netherlands. This chamber was the best in the world 
for providing thermal, vacuum and solar simulation 
with a dynamically quiet seismic mass boundary. As 
shown in Fig. I, SA3 consists of two wings. NASA 
placed one of the flight wings in the LSS on a stiff 
pedestal and rotated the wing 49 degrees to the solar 
beam. This rotation and the 2.9 m tall pedestal was 
necessary to position the wing for full illumination by 
the solar simulator with its 6 m beam. Fig. 2 shows the 
SA3 flight wing mounted on the pedestal. The 656 lb 
SA3 spans a width of 7.6 m and a height of 4,1 m 
above the top of the pedestal. 

The testing was conducted over four days using two 
shifts to cover 24 hours per day. After contamination 
bakeout and thermal equilibrium, the LSS simulated 20 
orbital cycles. During the orbit simulations, data 
acquisition systems measured 40 hours of signals from 
instrumentation. This paper describes the methodology 
used to screen the data to find any thermally induced 
vibrations of the SA3 wing. 

2 INSTRUMENTATION 

For this test, the instrumentation consisted of 
accelerometers, force transducers, thermocouples, a 
vacuum gage and solar sensors. Fig. 2 shows the 
location of the accelerometers with labels al 
through a7. All of the locations had three axes 
instrumented with the exception of a5. Only the X sa 
and Z sa axes were instrumented at a5. 


Fig. 1 HST with Solar Array 3 Installed 
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Fig. 2 Flight Solar Array 3 on Test Fixture 

Fig. 3 shows the locations of the four force transducers 
located at the top of the pedestal. They are labelled fl 1 
through fl 4. Each force transducer produced three 

orthogonal force signals oriented in the pedestal 
coordinate system. Two of the pedestal axes are 
labelled X ped and Z ped with X ped pointing away from the 
sun simulator. The 131° angle between the SA3 and 
pedestal coordinate systems is shown. Eight of the 
force transducer signals were col linear; therefore, the 
signals were summed before the charge amplifier. This 
reduced twelve force signals down to a total of eight. 



Fig. 4 shows the installation of the force transducers. 
They were sandwiched between an aluminum plate and 
the top of the stainless steel pedestal. A bolt through 
the center of the transducer was torqued while 
monitoring the force signal from the transducer. This 
technique allowed each bolt to be preloaded to 2000 lb. 



3 SIGNAL PROCESSING 

All of the primary and secondary bending and torsion 
modes of the SA3 exist between 1 and 10 Hz. 
Therefore, the data was band pass filtered in the 
following manner. An 8 th order digital filter removed 
background noise above 10 Hz from the measured 
data with parameters set to emulate an analog elliptic 
low pass filter. Then, a 4 th order digital filter removed 
transducer noise below 1 Hz with parameters set to 
emulate a Butterworth high pass filter. Application of 
these filters is described in [3]. 

The energy of the filtered signals y, may be described 
in relation to a mean square value, HK 2 , on p. 14 of [2]. 
This is simply the average of all the squared values in 
the time history as shown in Eq. 1. Table 1 contains 
filtered signal results of Eq. 1 for the two hour long 
13 th orbit simulation. Table 1 uses units of milli-pounds 
(mlb) and micro-g’s (pg). 
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In order to find short duration non- stationary events, 
Eq. 2 defines a short period average over time T to 
estimate the mean square values as a function of time. 
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Using the square roots of Eqs. 1 and 2, Eq. 3 
normalizes the root mean square time histories of the 
filtered signals. 
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Table 1 Orbit 13 Signal Root Mean Square (rms) 



4 SIGNAL SCREENING 

With 40 hours of data to process, a screening method 
for finding significant events of short duration was 
needed. The background noise level was fairly 
stationary with little variation from orbit to orbit 
simulation. This is evident in Fig. 5 showing the 
minimum, maximum and mean of all the orbit 
simulation power spectral densities (PSD) for the Y 
axis signal of force transducer f 14. 



Fig. 5 PSD of Force Transducer f!4 Y sa Axis Signal 


On p. 234 of [2], the reference suggests the use of short 
period mean square estimates to test for stationarity. In 
this analysis, time histories of normalized short period 
root mean square values were used in screening the 
data to uncover short duration events of non- 
stationarity. 

In the following sections, three different screens were 
developed. The first simply looked for any force 
transducer signal exceeding an arbitrary threshold over 
the noise floor. Next, the double screen added the 
requirement that two force signals and two 
accelerometers must all exceed a threshold at the same 


moment in time. This is reasonable, since any primary 
and secondary mode wing vibrations would result in 
multiple responses. Finally, the triple screen used the 
boundary accelerations measured at the top of the 
pedestal to analytically predict the force transducer 
signals resulting from chamber vibrations. Whenever 
the predicted force signals exceeded a threshold during 
an event found in the second screen, the event was 
removed from consideration. 

4.1 First Screen 

After computing , these results were screened to 
find events that might cause HST to have a loss of lock 
event. For this screen, only the force transducer 
signals were examined. The screen flagged each time 
when any one of the signals had an upward crossing of 
thresholds. This screen can be described with Eqs. 4-7. 
Eq. 4 generates a threshold flag time history of ones 
when T' exceeds s and zeros otherwise. Eq. 5 sums all 
the force transducer threshold flag time histories to 
produce a time history that provides a count of all the 


force signals that exceeded Eq, 6 generates a screen 
logic time history of ones for the times when any one 
of the force signals ¥ exceeds r. Finally, Eq. 7 

produces an up tick time history of only the rising 
transitions of p\. 
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The up tick p e i time histories were summed for each of 
the orbit simulations. Table 2 lists a summary of the 
results of this screen. For a r of 4, a total of 25 events 
were found. These events were identified during the 
test using another methodology and were analyzed by 
HST pointing and control engineers. Their pointing 
and control system simulation results predicted no 
LOL for HST. This was very good news, and as a 
result the LOL disturbance requirement for SA3 was 
verified. In addition, the test objective was met. 

In Table 2, orbit 13 was observed to be more active 
than the other orbits. Th is was also noticed during the 




15 / 


test. At the time, the engineers hypothesized that 
during an earlier orbit the activati on of a heater on the 
aluminum interface plate, shown on Fig. 4, was the 
cause. Based on this hypothesis the test director 
disabled the heater. 

Table 2 Maximum Force rms Events 
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Fig. 6 contains time history plots of the maximum 
normalized rms signals grouped by SA3 accelerometers 
in the top plot, force transducers in the middle and 
boundary accelerometers in the bottom plot. The eight 
events counted in Table 2 for orbit 13 with a z of 4 can 
be seen clearly in the middle plot of Fig. 6. Notice that 
the major peaks on the top plot do not line up with the 
major peaks of the middle plot. Since the 
accelerometer and force transducer peaks did not 
correlate, they cannot be the result of SA3 vibrations. 

4.2 Double Screen 

Although the first screen identified no events that 
would cause LOL, any thermally induced vibrations 
were still interesting for their jitter effects. In a further 
attempt to find SA3 vibrations, this double screen was 
devised to find events with at least two forc e 
transducer signals and at least two SA3 accelerometer 
signals greater than r. Eqs. 4-9 define the screen. Eq. 8 
sums all the SA3 accelerometer signals from al to a5 
that have a q* greater than r. Eq. 9 defines the screen 

logic time history. 
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Fig. 6 Orbit 13 Normalized Signal rms 
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Table 3 summarizes the results of the double screen. 
The number of events with a z of 4 was reduced from 
25 to only 2 events. All of these events occurred 
during the 5 th orbit simulation. 

Table 3 Double Screen Events 




E p??<, 2 

) with T= Is | 

orbit 

2.0 

2.5 

3.0 

3.5 

4.0 

1 

42 

3 




2 

64 

7 




3 

41 

3 




4 

68 

2 




5 

70 

13 

4 

3 

2 

7 

87 

10 

3 



8 

82 

6 




9 

94 

10 

4 



10 

45 

7 

1 



11 

74 

4 




12 

62 

5 




13 

92 

6 




14 

77 

4 




15 

79 

6 




16 

74 

4 




17 

82 

6 




18 

86 

5 




£ 

1219 

101 

12 

3 

2 





Fig. 7 plots the maximum % time histories for 5 th orbit 

simulation. Again, the top plot is the maximum of all 
the SA3 accelerometers. The middle is the force 
transducers, and the bottom plot is the seismic mass. 
The two events with a r of 4 are clearly evident on all 
three plots at 3500 and 4000 seconds. Since the seismic 
mass was moving, the SA3 responses were the result of 
base motion. 


For the final triple screen, the predicted force 
transducer responses were used in Eq. 1 1 to sum the 
number of W pf greater than z p . Eq. 12 defines the 

logic time history. It is the same as the double screen 
with the added requirement on the predicted forces jy . 
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Fig. 7 Orbit 5 Normalized Signal rms 


This assumption was verified by simulating impulse 
responses of the SA3. In the simulations, responses 
were observed on the accelerometers and the force 
transducers, but they were not observed on the top of 
the pedestal. 


4.3 Triple Screen 

In order to continue the search for thermally induced 
vibrations of the SA3, this third and final screen 
eliminates boundary excitation events by using a test 
data derived linear model to predict SA3 responses, A 
system identification of SA3 using the test data was 
performed and resulted in the A, B, C and D matrices 
contained in Eq. 10 for computing linear predictions of 
the SA3 responses y . The methodology for 
performing the system identification is detailed in 
section 5. 
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Table 4 lists the results from the triple screen. The 
number of events with a z of 2.5 was significantly 
reduced from 101 to a total of 4. On closer inspection 
of the time histories for these 4 events, they still appear 
to be caused by boundary motion. Therefore, no 
thermally induced vibrations of the SA3 were detected 
at z of 2.5. 


Table 4 Triple Screen Events 
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5 SYSTEM IDENTIFICATION 

The triple screen above used predicted SA3 responses 
resulting from base excitation. In order to perform 
accurate predictions, a system identification was 
performed to estimate the state space matrices in 
Eq. 10. The elements of the matrices were cast into a 
modal parameter grey model form to be used in a 
parameter estimation algorithm detailed in [4]. 

Using the accelerometer signals from location a6 as an 
input source and all the other signals as output, the 
modal parameters were extracted from the test data. 






Table 5 lists natural frequency and percent critical 
damping coefficient results from the system 
identification for both ambient and vacuum conditions. 
The right two columns contain the analytical frequency 
values from the finite element model of the test 
configured SA3 and critical damping coefficients from 
SA3 modal survey data for comparison. 

Table 5 Modal Parameters 


Since no interior forces f a are applied the second row 
of Eq. 16 may be written as Eq. 17. If a thermally 
induced vibration event occurred, this assumption 
would be false momentarily; however, the averaging 
described in the next section will still make this a 
reasonable assumption. 

l' = -2 £a>J-a> 2 n Z-m ih x h (17) 



Ambient 

Vacuum 

Predicted | 


fn 

c 

fn 

c 

fn 

G 

Mode 

(Hz) 

(%) 

(Hz) 

(%) 

(Hz) 

(%) 

1 

1.2 

3.1 

1.3 

4.7 

1.1 

4.7 

2 

1.5 

4.4 

1.7 

4.6 

1.4 

3.9 

3 

- 

- 

1.8 

3.0 

1.8 

2.3 

4 

5.3 

0.2 

5.5 

0.4 

5.2 

1.7 

5 

6.2 

0.9 

6.4 

1.4 

5.9 

1.4 

6 

6.3 

0.9 

6.6 

0.6 

6.5 

1.7 

7 

7.4 

1.0 

7.8 

0.7 

7.5 

1.6 

8 

12.9 

0.7 

13.7 

1.9 

12.9 

2.4 

9 

13.9 

2.7 

14.6 

0.8 

13.4 

2.9 


5.1 Modal Parameter Grey Model 


For the SA3, the boundary may be considered 
kinematic. Therefore, k bb is zero, and the first row of 
Eq. 16 may be written as Eq. 18. Eq. 19 is formed by 
taking the second derivative of the second row of 
Eq. 15. 

f h =m h J + m hh x h (18) 

+ (19) 

Finally, Eqs. 20 and 21 define the grey model state 
space matrices from Eqs. 17-19. 


Eq. 13 represents the physical motion for the SA3. The 
degrees of freedom can be partitioned as shown in 
Eq. 14, where the Xb is the boundary and X a are the 
interior degrees of freedom. 

Mx + Cx + Kx = F (13) 
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Substitute Eq. 15 into Eq. 14 and postmultiply by the 
transpose of the transformation. This Craig-Bampton 
transformation is described in [5], 
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5.2 Spectral Analysis 
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To estimate the modal parameters in Eq. 21, a spectral 
analysis of the data was performed. The objective was 
to compute frequency response functions using the 
boundary accelerations as a reference. 


Assuming diagonal damping and no boundary damping 
coupling, Eq. 16 shows the results of the 
transformation. 
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The power spectral density functions and cross spectral 
density functions in Eqs. 22 and 23 may be calculated 
from the cross-correlation function in Eq. 24, as on 
p. 435 in [6]. 
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In practice, the cross-correlation function in Eq. 24 is 
not calculated, and time cannot extend from-oo to+oc . 
An approximation of the spectral Eqs. 22 and 23 is 
made by averaging methods detailed in [3] and [7]. 

Now the frequency response functions may be 
computed with Eq. 25, as in p. 162 of [2]. For this 
analysis, the input degrees of freedom, denoted by the 
subscript u, were the signals from accelerometer a6 
shown in Fig. 2. Since only the translational 
accelerations were measured, S ^ is a three by three 

matrix at each frequency, CO . 

H )V (co) = S;l{co)S yv (c0) (25) 

Eq. 26 computed the multi-input coherence C y of the 
frequency response functions. The operator on the third 
term of the numerator in this equation takes transpose 
of the complex conjugate of H }V . Eq. 27 computed the 
phase of the frequency response functions (FRF) by 
taking the imaginary part of its natural logarithm and 
converting to degrees. Fig. 8 shows a FRF plot of the 
force transducer fl4 Y sa axis signal referenced to the 
Z,* axis of accelerometer a6. 
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5.3 Modal Parameter Estimation 

This section will describe how the parameters were 
computed to form the fit shown in Fig. 8. Using the 
FRF’s computed by Eq. 25, impulse response functions 
in Eq. 28 were calculated by using the real portion of 
their inverse Fourier transformation. Then, the impulse 
response functions were truncated and trained together 
in Eq. 29. For the input signals, three impulses are 

defined in Eq. 30 where <^is the Dirac delta function. 



Fig. 8 Orbit 13 Force Signal FRF and Coherence 
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Fig. 9 contains a plot of the impulse response functions 
computed for the three FRF’s of the Y sa axis force 
transducer signal from location fl 4 referenced to the 
three accelerometer signals from location a6. 



The state matrices in Eq. 21 are a function of the 
parameters contained in the parameter vector defined in 
Eq. 3 1 . Here the elements in matrices <f>, m bb , m K and 
T ba are reordered into single dimensioned vectors. 
Initial values for £ co„ and <j> are required to start the 
parameter estimation process. 
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To obtain starting values of Q and (On, a subspace 
method estimated a black box of the A, B, C and D 





matrices with u and y ( as objective inputs and outputs, 
as in [4]. Then the complex eigenvalues, A, of A were 
computed. The natural frequencies, cd„ , were computed 
from the magnitude of A, and the critical damping 
coefficients, were computed by taking the negative 
cosine of the imaginary portion of the logarithm of A in 
Eq. 32. 

co n =\A\ Q = -cos(Im(InA)) (32) 

From these values, only reasonable results were 
retained. Then, Eq. 33 computed an initial value of <f> 
by using the imaginary part of the peaks where j was 
selected to give the largest response at the natural 
frequency of The remaining parameters in 0 may 
be started at zero. 

4=I*n (33) 

The state space matrices in Eq. 21 were estimated with 
a parameter estimation routine that minimizes the error 
between the predicted time history y and the impulse 
response functions y , detailed in [4]. Results of the 
parameter estimation are shown in Fig. 8 and in 
Table 5, and they are considered to be very good. 

6 CONCLUSIONS 

The disturbance verification test objective was to 
measure thermally induced vibrations of SA3 in a test 
simulated space environment. The events of interest 
were transient vibrational responses. It was mandatory 
to measure these vibrations if they had been large 
enough to cause Loss of Lock (LOL) on the HST guide 
stars. [1] 

The objective of the test was met. Events were 
identified in the data. However, based on HST pointing 
control system simulations, no measured disturbances 
were large enough to cause the HST to experience 
LOL. Furthermore, many of the events identified by 
the first screen were not attributed to thermally induced 
vibrations of the SA3, since the accelerometer and 
force transducer signals did not correlate. 

The double screen added the requirement that two 
force signals and two accelerometers must all exceed a 
threshold. This requirement was an attempt to 
eliminate noise on a single transducer signal from 
triggering and event. This screen found many events 
that were easily identified as boundary motion, since 
boundary motion will generate SA3 responses on both 
the force transducers and accelerometers. 


Finally, the triple screen used the boundary 
accelerations to predict the force transducer signals. 
Whenever the predicted signals exceeded a threshold, 
the events were removed. The results of this final 
screen essentially eliminated all the events that were 
above the background levels. In other words, no 
thermally induced vibrations of the SA3 were detected 
above the background noise. 

At this point in the analysis of the test data, no 
significant Solar Array 3 induced disturbances have 
been found in the test data that would impact the 
performance of the Hubble Space Telescope. 
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